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ABSTRACT 

We investigate under which conditions cold, fan-shaped winds can be steadily 
launched from thin (Keplerian) accretion discs. Such winds are magneto-centrifugal 
winds launched from a thin annulus in the disc, along open magnetic field lines that 
fan out above the disc. In principle, such winds could be found in two situations: (1) 
at the interface between an inner Jet Emitting Disc, which is itself powering magneto- 
centrifugally driven winds, and an outer standard accretion disc; (2) at the interface 
between an inner closed stellar magnetosphere and the outer standard accretion disc. 
We refer to Terminal or T-winds to the former kind and to Magnetospheric or M- winds 
to the latter. 

The full set of resistive and viscous steady state MHD equations are analyzed 
for the disc (the annulus), which allow us to derive general expressions valid for both 
configurations. We find that, under the framework of our analysis, the only source 
of energy able to power any kind of fan-shaped winds is the viscous transport of 
rotational energy coming below the inner radii. Using standard local a prescriptions 
for the anomalous (turbulent) transport of angular momentum and magnetic fields 
in the disc, we derive the strength of the transport coefficients that are needed to 
steadily sustain the global configuration. It turns out that, in order for these winds 
to be dynamically relevant and explain observed jets, the disc coefficients must be 
far much larger than values expected from current knowledge of turbulence occurring 
inside proto-stellar discs. 

Either the current view on MHD turbulence must be deeply reconsidered or 
steady-state fan-shaped winds are never realized in Nature. The latter hypothesis 
seems to be consistent with current numerical simulations. 

Key words: Accretion, accretion discs - Magnetohydrodynamics (MHD) - stars: 
formation - stars: mass loss - stars: pre-main sequence - ISM: jets and outflows 



1 INTRODUCTION 

1.1 Mass loss from young stars 

Young Stellar Objects (hereafter YSOs) show evidences of 
fast (several hundred km/s) and collimated jets that carry 
a sizable fraction of the released accretion power (e.g. |Bally| 
et al. 2007 Ray et aL||2007| ). It has been soon recognized 
that large scale magnetic fields play a major role in driv- 
ing these jets. Basically, this is the only process that has 
been proven to be able to simultaneously provide a very 
efficient acceleration along with a self-confinement of the 
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ejected plasma (see Cabrit 2007 for a review). These large 
scale magnetic fields can be anchored onto the star (giv- 
ing birth to a "stellar wind", e.g. [Hartmann & Mac Gregor| 



19§2}|Sauty fc Tsinganos|p94l [Matt fc Pudritz||2005|), the 



accretion disc ("disc wind", e.g. Blandford & Payne 1982 
Pudritz fc Normanl [19861 |Wardle fc Konigl||1993| |Ferreira 



Shu et al. 



fc Pelletier||1993a[) or linked to both ("X- winds' 
1994 a||Cai et al.|2008| and " Reconnection X- winds" , |Ferreira 
|et al.|2000| >. 



Each model corresponds actually to a different scenario 
for the origin of the magnetic field. On one extreme, disc 
winds rely on the existence of a large scale magnetic field 
threading a significant portion of the disc (up to ~ 100 AU 
for instance in earlier works). On the other extreme, stel- 
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lar winds and X-winds assume the presence of such large 
scale fields only on the star. Since direct observations of 
magnetic fields in discs are highly difficult (see for instance 
Donati et al^|2005| |, some conjectures must be made. How- 
ever, these wind models are described by the same set of 
magnetohydrodynamic (hereafter MHD) equations and can 
be easily confronted to optical and near-IR observations of 
YSO jets (Ferreira et al. 20061. For instance, if the trans- 



verse velocity shifts observed in some of these jets are indeed 
due to jet rotation, then stellar winds and X-winds (in fact 
any wind coming from the innermost regions) are ruled out 
as the main jet component: they simply do not carry enough 



angular momentum to explain observations ( Anderson et al. 
2003 Ferreira et al.|2006 K But on the other hand, these op- 



tical observations imply also that the last relevant field lines 
in disc winds cannot be anchored farther than a few AU (for 
the atomic component: molecular jets could be anchored up 
to ten AU, |Panoglou et al.||2012| >. Thus, whether YSO jets 
are rotating or not, observational evidences suggest that jets 
are arising from regions smaller than or of the order of sev- 
eral AU, while discs could be as large as several hundreds of 
AU. 



1.2 What are fan-shaped winds? 

If observed YSO jets are indeed mainly fed by disc winds, 
then the mass, angular momentum and energy extraction 
due to the winds have a profound impact on the dynamics 
(and possibly observational appearance) of the underlying 
Jet Emitting Disc (hereafter JED, Combet fc Ferreira|2008 



Combet et al.|2010[). It has been analytically shown (|Ferreira| 
fc Pelletier||1995||Ferreira|1997||Casse fc Ferreira|2000| > and 



numerically confirmed (Tzeferacos et al. 2009[) that the large 



scale magnetic field must be close to equipartition in order 
to magnetically launch jets in a steady fashion. A natural 
transition from an inner JED to an outer standard accretion 
disc will therefore occur if, beyond a radius r j , the local disc 
magnetization suddenly decreases. If this drop is done over 
a thin radial extent, namely Ar ~ h <C r where h is the 
local disc vertical scale height, then the poloidal magnetic 
field tends to take a fan-like geometry above the disc (see 
Fig.[l^). This region of highly bent, open magnetic field lines 
could also be the locus of some ejection, allowing thereby to 
act as a transition region between the outer non ejecting 
disc to the inner JED. These MHD winds will be hereafter 
referred to as " Terminal- winds" or T- winds. 

A different sort of fan-shaped MHD wind could also be 
found at the interface between an outer standard accretion 
disc (with no or negligible large scale magnetic field) and the 
stellar magnetosphere. For simplicity, let us consider only 
the case of a stellar magnetic moment aligned with the ro- 
tation axis of the star and the disc. The disc is truncated at 
a radius R t where the stellar magnetic field is near equipar- 



tition with the disc thermal pressure 1 


Romanova et al.|2002| 


Bessolaz et al. 2008 . Zanni & Ferreira 


2009||Romanova et al.| 


2011 


2012 


1 . Such an interaction leads to the opening of stel- 



lar field lines so that some magnetic flux diffuses in the disc. 
Because of the sharp decrease of the (initially) stellar dipole 
field, the zone where this field is dynamically important is 
very thin, a typical length being Ar of order of a few scale 
heights. Here again, such a geometry appears to be a favor- 
able place for magnetically launching winds (see Fig. ^p)- 



These " Magnetospheric- winds" or M- winds are anchored in 
a radial zone that could act as a buffer zone between the 
outer non-ejecting disc and the inner magnetospheric accre- 
tion curtains. 

In the case of M- winds, three different situations can 
take place depending on the location of the disc truncation 
radius R t with respect to the co-rotation radius R co (defined 
as the radius where the Keplerian Qk and stellar f2* angu- 
lar velocities are equal). When Rt < R co accretion columns 
will be formed whereas a " propeller regime" is settled when 
Rt > Rco (|Romanova et al.|2003||2005||Ustyugova et al.|2006 



and references therein) . The limiting case where Rt — R co is 
known in the literature as "X-winds" and its axisymmetric, 
steady state ideal MHD wind structure has been explored 



in several papers (Shu et aL]|1994a|b| Najita & Shu 1994 



Shu et al.lpkl |Shang et al.f|1998| 2002| |Cai et al.||2008 l 



The computed 2D flow structure is a fan-shaped cold wind 
expanding from a point source with an initially super slow- 
magnetosonic speed that is matched to an asymptotic jet 
cylindrical solution. In these wind calculations the under- 
lying disc is a mere boundary condition and the impact of 
the winds on the underlying disc dynamics has not been 
investigated yet. 

The goal of this paper is to analyze at which condi- 
tions fan-shaped winds, namely T- winds and M-winds, can 
be launched from a near Keplerian accretion disc. By fan- 
shaped wind, we mean any steady-state axisymmetric MHD 
wind that originates from a region of extent Ar ~ ft in a 
Keplerian accretion disc and with the tendency to fill in all 
space (unless confined by an outer pressure). The questions 
we plan to address are the following: how much mass and 
power can be carried away by these winds? Is a steady state 
possible? What are the constraints imposed on the underly- 
ing MHD turbulence in the disc? 

In Section 2, we write down and analyze all steady- 
state MHD equations describing the disc and its winds. Note 
that the disc is actually a thin annulus, from r^ to r e = 
n + Ari with Ar; ~ h. Our results will be general, valid for 
any kind of fan-shaped wind. Then, in Section 3 we apply 
them to each specific case, namely T-winds and M-winds, 
and address the questions raised above. The caveats and 
consequences of our analysis will be discussed in Section 4. 



2 GENERAL ASPECTS OF DISCS DRIVING 
COLD, FAN-SHAPED WINDS 

In this Section, we write down all equations and put them 
into a form that is general and valid for the two envisioned 
fan-shaped wind scenarios, T-winds and M-winds. Next Sec- 
tions will be devoted to each of these scenarios by applying 
our general results to each specific case. 



2.1 Governing MHD equations for the disc 

To describe the disc, the following set of non-relativitistic, 
axisymmetric, steady-state single fluid MHD equations must 
be analyzed 

V-pu = (1) 
pu ■ Vu = -pV$ G - VP + J x B + V ■ T (2) 

VmJcP = Up X Bp (3) 
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Magneto-centrifugally driven Wind 




Jet emitting disc (JED) Standard disc (SAD) Stellar magnetosphere Standard disc (SAD) 



Figure 1. The two configurations of fan-shaped winds. Left: Terminal or T-winds launched from r 4 = rj, the outer radius of a JED, 
to r e = r; + Ari the inner radius of a standard accretion disc (SAD). The T-wind leans against the inner magneto-centrifugally driven 
wind and fans out above the standard disc. Right: Magnetospheric or M-wind at the interface between an outer standard disc and the 
inner magnetosphere of a central star, with Ti ~ r^f, where r^f is the base of the funnel flow |Bcssolaz c t al,|2008} . In this configuration, 
three cases can be realized according to the strength of the stellar magnetic dipole as measured by the ratio of the truncation radius 
Rt < Ti to the co-rotation radius R co '. M- winds with n < R C o, M- winds with ri > R co and X-winds, which are M- winds with n ~ Rco- 



V -(^VrB ) = V -^{B^Up-BpQr) (4) 

where J = V x B//i is the electric current density, 4>g = 
—GM/(r 2 + z 2 ) 1 ^ 2 is the gravitational potential of the cen- 
tral object and other quantities have their usual meaning. 
Here, T is the turbulent stress tensor resulting from a self- 
sustained MHD turbulence within the magnetized disc. Fol- 



lowing Shakura & Sunyaev (19731, the only relevant com 



ponent used here is T r< p = pv v r \ ~q~\ ■ ur mean field ap- 
proach assumes that the disc turbulence provides anomalous 
transport coefficients such as a turbulent viscosity v v and 
magnetic diffusivity v m (rj m — [i v m is the resistivity). This 
diffusivity, appearing in the toroidal component of Ohm's 
law Q , allows the diffusion of mass across the poloidal mag- 
netic field Bp and thereby the possibility of a steady state. 
A diffusion process must also be present with respect to the 
toroidal motion, otherwise the field would be overly sheared 
with a tremendous toroidal field component B^. In the in- 
duction equation Q, some turbulent magnetic diffusivity 
vL has to be introduced, so that the dimensionless number 



X™, = 



(5) 



measures a possible anisotropy of the turbulent diffusion 
processes within the disc ( |Casse fc Ferreira| 2000). The am- 
plitude of the turbulent transport coefficients are measured 
using ct coefficients, namely ct v — v v /Q.kh, ot m — v m /Q.Kh 
and a' m = v' m /Q.Kh, so that the effective magnetic Prandtl 
number writes 



On 



(6) 



These a coefficients are expected to decrease vertically on a 
disc scale height, so that the flow (the wind) becomes invis- 
cid and described by ideal MHD equations. It is noteworthy 
that these a coefficients can also vary radially. Indeed, their 
magnitude depend on, e.g. the field strength and the degree 
of ionization, both effects that vary throughout the disc. For 



instance, one might have the viscosity parameter ct v going 
from a very low value at the disc outer edge (cold disc and 
low magnetic fields) to a much larger value (possibly reach- 
ing unity) at the disc inner edge. On the other hand, while 
the degree of ionization has a critical impact on the develop- 
ment of MHD turbulence (see eg. |Lesur fc Longaretti|2007| 
Bai fc Stone||20l"T and references therein), its radial depen- 
dency is unknown and depends on the assumptions made 
about the disc density structure (see for instance |Fromang| 



et al. 2002 Combet et al. 20101. To circumvent these un- 
certainties, we will hereafter assume that V m undergoes no 
significant change within the portion of the disc considered 
here. This is clearly a limitation of our model. 

Finally, in order to close the above set of equations, 
we will assume a perfect gas law in the disc, namely P — 
pC 2 , with a vertical isothermal approximation. This latter 
assumption is not critical for our purpose as we deal with 
cold ejection from near Keplerian accretion discs. 

The dynamical importance of the large scale magnetic 
field is measured at the equatorial plane by the disc magne- 
tization 



where Ptot is the total (gas + radiation) pressure in the disc 
( Ferreira fc Pelletier|1995| . In a disc with a negligible radi- 
ation pressure as considered here, fj, = 2/(3 where /3 is the 
usual beta plasma parameter. Maintaining a disc with a Ke- 
plerian profile despite the presence of large scale magnetic 
fields puts an upper limit on the magnetization: p cannot be 
larger than unity (Ferreira & Pelleticr 1995 Ferreira||1997 



Shu et al. 20081. Under these circumstances, the magnetic 



pinching force is at most comparable to the gravitational 
component and the disc aspect ratio can still be approxi- 
mated by 



P n 2 K r* 



C s 



< 1 



(8) 
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namely C 3 ~ Qnh, as in a standard accretion disc. The 
actual scale height of the disc will be smaller than, but of 
the order of h, which is enough for our purpose here (see 
Appendix [B] for more details). 



2.2 Cold winds and magnetic field configuration 

Fan-shaped winds are assumed to be cold, that is magnet- 
ically driven. This implies that the magnetic force must be 
the dominant one above the disc and that enthalpy can be 
neglected. A successful magnetic ejection then requires some 
conditions at the disc surface that can be translated into con- 
ditions on the magnetic field geometry, namely on the radial 
and toroidal B^ magnetic field components at the disc 
surface. 

The energetic condition to launch cold jets from near 
Keplerian discs has been first derived by Blandfor d fe Payne| 
(1982). It states that the poloidal magnetic field lines at 



the disc surface must be inclined by more than 30 deg with 
respect to the local vertical. This translates into 



P = 



(9) 



larger than, but of the order of unity, where B z is assumed 
to remain of the order of the field at the midplane (see Ap- 
pendix [A!. One simple way to look at the ejection process 
is to project the magnetic force along a magnetic poloidal 
surface, which gives 

F„ = (J x B) • e = ^V,|7 



(J x B) • e|. 



B^ 
2nr 



Vi,/ 



(10) 



where B p = B p e^, Vii = eii • V and I = 2-ktB^I '(j, < is 



the total current flowing within this magnetic surface ( Fer- 
|reira|19 97). These expressions show that plasma acceleration 
along any given magnetic surface can be seen as a current 
leakage through this surface, providing simultaneously an 
azimuthal and a poloidal acceleration. In other words, while 
the magnetic field is both braking down (f^ < 0) and verti- 
cally pinching (F\\ ~ F z < 0) the disc, the sign of these two 
forces must change at the surface so that ejection becomes 



possible ( Ferreira fe Pelletierj 1 993b 1995|. This can be nat- 
urally achieved but imposes a constraint on the magnetic 
shear so that 



f_ 

B z 



(11) 



This important condition is derived from the MHD induc- 
tion equation Q applied in a Keplerian disc where the dom- 
inant contribution is B r dQ/dr (because of the magnetic ge- 
ometry considered here, see Appendix [c] for details). 



2.3 The disc angular momentum transport 

The disc angular momentum equation (^-component of 
Eq. [2| writes 

r r* or 

where F$ = J z B r — J r B z . Integrating this equation over the 
disc thickness gives 



_ 3 m 

dr 



(12) 



— ~ -2 h ——T,v v Q,Kr (13) 

vkt p a 2r 2 or 

where both O ~ Qk and h roughly constant over the extent 
An have been assumed. Here, £ stands for the local disc 
surface density, defined as £ = J pdz. Defining the accre- 
tion sonic Mach number as m 3 — M a /2iiT,^lKrh — u /C s , 
where u is the average accretion velocity, one gets 



+ m s 



2qp + (a v e 



(14) 



where m s na3 and m" lsc are respectively the magnetic and 
viscous contributions and where 

_ dW^vy^) _ 3 riln7l e d\nM a 

dlnr - 2 6 d\nr +6 dlnr ( °> 

Here, the vertically averaged Reynolds number of the accret- 
ing flow writes TZ e — ru a /v v — m s /a v e = £(1 + A), where 
A = m™ as /m" ,f " : is the ratio of the magnetic to the vis- 
cous torque. In a standard accretion disc (hereafter SAD), 
the magnetic torque is negligible and lZ e = (" = 3/2. In a 



JED, where the accretion rate varies as M a oc (Ferreira 



fc Pelletier|1993a| |, one gets £ = 3/2 + 3£. Note also that, in 
a steady-state picture, the viscous torque is always braking 
down the disc plasma unless lZ e increases with the radius. 
In usual accretion disc models, TZ e is however a constant. 



2.4 Properties of cold fan-shaped winds 

Above the disc, turbulence is assumed to decay and ideal 
MHD applies. Steady-state, axisymmetric cold MHD jets 
can be described with 4 invariants: magnetic flux to mass 
flux ratio rj, surface angular velocity Qb, total specific an- 
gular momentum L — Q.bTa (where ta is the Alfven cylin- 
drical radius) and specific energy E. Our objective in this 
section is to relate these invariants to the disc physics, al- 
lowing henceforth to link cold fan-shaped wind asymptotic 
properties to disc properties. 

Mass and magnetic flux conservation provide 



flop 



B,; 



(16) 



where ri(<&) is a constant along a magnetic surface defined 
by a constant $ = J 2ixrB z dr. It measures the ratio of the 
mass flux to the magnetic flux carried in by the wind, namely 
M w = J 2-Krpu z dr = J — where the integration is made 
over the total magnetic flux $x leaving the region of extent 
An. 

The induction equation Q becomes in ideal MHD 



P-oPr 



(17) 



where f2* is the rotation rate of a magnetic surface, assumed 
here to be the Keplerian value at r;. Because of plasma in- 
ertia, magnetic field and plasma have not the same rotation 
rate, which induces a toroidal field component. Plasma ve- 
locity can be written u = (^//^ p)B + £!*re^ and is not 
parallel to the total magnetic field. In the reference frame 
rotating at SI* however, they do become parallel. 

The only explicit calculations of stationary fan-shaped 
MHD flows that can be found in the literature are those re- 



lated to the X- wind ( jNajita fc Shu||1994| |Shu et al.||1995 
Cai et al.||2008[ ). To ease comparison with these studies, we 
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will use hereafter the same notations, with Rx = f"i. Defin- 
ing fiducial quantities Mx = M w , px = Mx /2irrx&x , 
Vx = VxRx, B x = (no flxMx/2nR x ) 1/2 _ and Q x = 



Shu et al. 



(1994a I, the dimen- 



(GM/R 3 X ) 1/2 accordingly to 
sionless velocity and magnetic field are related by B 
/?($)pu where 



PoPxVx 
r/B x 



-i Bx 



(18) 



In this expressio n, kbp = r)Q,K,iTi/ 'B Z:i is the Blandford 



& Payne ( 1982 1 mass loading parameter. Using M u 



J dM w = fdWpxVx/PBx one gets 



2-kBxRx = 



d<t> 



(19) 



10 

where j3 is some average value in the wincQ Since this mag- 
netic flux is arising from a tiny disc region around a radius 
Ti, we can estimate it as $x — 2,rB z ^riAri, where B Zt i is 
the vertical magnetic field. Its value is then related to the 
fiducial field Bx by B z ^ ~ fiBxTij Ari. This is the cor- 
nerstone of fan-shaped winds: a spherical dilution of the 
magnetic field, giving an average mass loading parameter 
Rbp = fi~ 2 Arj/Vj. Using the definition of the disc magneti- 
zation at Ti brings 



Pi ~ P 



2 ( 2Af„ 



ArJ 



f 



/ A) 7 



(20) 



which is general, valid for any fan-shaped wind configura- 
tion. In this expression, we introduced / = as the ratio 
of the total (two-sided) mass loss in winds to the incoming 
mass rate in the disc at r e . 

The angular momentum conservation in the jet writes 
L($) = Q,,r 2 A = Q,r 2 ~ rB^/r/ where r A is the Alfven radius 
and L(&) the total specific angular momentum carried along 
a magnetic surface. Normalizing it to Q.K,ir 2 and evaluating 
it at the disc surface gives 



&K,ir 2 



q 

KBP 



An 



(21) 



which is valid for any fan-shaped wind configuration. Note 



that J = Asp, the Blandford & Payne ( 1982 1 magnetic lever 
arm parameter. 

The projection of the momentum conservation equation 
along a magnetic surface provides the Bernoulli invariant 



rBi, 



(22) 



which is the total specific energy carried by each mag- 
netic surface. However, only a fraction £(<£>) is actually 
available for accelerating the plasma, the part Q 2 r 2 A (the 
energy of the magnetic rotator) remaining unused. If, at 
infinity, the magnetic field structure keeps almost no en- 
ergy, then the flow reaches its maximum allowable speed 
= n Kl n(2J($)-3) 1/2 . 
The total power carried away by the cold wind through 
the disc surface is 



1 If we define $ 
fy)- 1 / 3 used by 
by 9/8. 



and take the scaling /?(*) = |/3(1 - 
1 2008}, the rhs of Eq.ll9| is multiplied 



Cai et al 



2P wi nd = 2 E{$)pu p ■ dS = 2Pmhd + 2P kin 



(23) 



The kinetic power carried by the outflowing plasma at the 
disc surface, namely 



Pkin = 



U -r \ ,„ GMM a 

i + * G P u,dS,-/ — 



(24) 



is initially negative: the wind power is mostly stored as mag- 
netic power that will be eventually transferred to the plasma. 
The dominant contribution in P W md is the MHD Poynting 
flux S M hd leaving the disc 



MHD 



<>,T 



m 



Po 

GMM ai 2gp 
4r; m s 



-B p • dS 



B a B z 



Po 



Ar, 



2-KTi O^Ar, 



(25) 



Inserting Eq. ( 20 I into this expression provides a wind power 
GMM ae r , 2qp 



2-Fwind 



2r. 
GMM a 



An 



((1-/) 
2n \ An ) 



') 



(26) 



Defining the average maximum jet velocity 

Vui aS Pwind 

\M w v%, one obtains 



v w = VL K iTi ( 2g/3 



An 



1/2 



(27) 



This expression is nothing more than v w = ^KiTi(2J— 3) 1,/2 , 
with an average magnetic lever arm J = l + qj3 2 ri/ Art. This 
result is again general and valid for any cold fan-shaped wind 
arising from a near-Keplerian accretion disc. 



2.5 The global energy budget of the disc 

The local energy conservation law for a steady-state MHD 
flow is given by the following expression 



pu p 



+ $c 



+ S rad - u • T = 



(28) 



where $g, Smhd, S ra d, T are respectively the gravitational 
potential of the central star, the MHD Poynting flux, the 
radiation flux and the turbulent "viscous" stress tensor (u 
stands for the velocity of the fluid and p for the plasma 



density). The last term in Eq.(28l mimics the energy trans- 
port by turbulence and vanishes at the disc surface. The 
plasma enthalpy may be neglected here. Using the Green- 
Ostrogradky theorem on the volume defined by the disc (ac- 
tually an annulus of extent Ari and height h) allows to de- 
rive the global energy budget, namely 



+ P vi 



(29) 



In this equation, the total power carried away by the cold 
winds 2P W i na (sum of kinetic and magnetic energy fluxes) 
and the disc luminosity 2P ra( ; are fed by the power available 
in the disc between n and r e . This power has two possible 
origins: the released accretion power P acc and the power P v i 3 
brought in by viscosity. Any other energy flux, like e.g. heat 
deposition by irradiation, is assumed negligible here. 
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2. 5. 1 The accretion power 



2.5.3 The disc luminosity 



By definition, the accretion power is the mechanical power 
released by the plasma accreting from the external radius of 
the launching zone r e to the inner radius n, namely 



+h 



+ <3?g I 2nrpu r dz 



GMMa 

2r 



(30) 



where [A]^ — A(r e ) — A(ri). Because of the wind, the disc 
accretion rate M a (r) = — J , 2nrpu r dz is a function of the 
radius. Defining M ai — M a (ri) and M ae = M a (r e ), mass 
conservation writes M a i = M ae — 2M W = M ae (l — /), leading 
to 



GMM ae (An 



2r, 



(31) 



It can be readily seen that there is no mechanical power 
available if / > An/n- there is not enough radial contrast 
to feed large amounts of mass with mechanical power. 



2.5.2 The "viscous" power 

In a turbulent accretion disc, the inward flux of mass co- 
exists with outwardly directed fluxes of energy and angular 
momentum. At each radius, energy and angular momentum 
are therefore being deposited from the inner radii, which 
translates here into a net power 



p 

i vis — 



(u • T) • dS 

i GMM, 



2-KrQrT r &dz 



.! GMM a ,e 



where 
23 = (1 



/) 



7^ 

Tie 



(32) 



(33) 



Since the only available sources of energy are P acc and P v i s , 
the difference factor T> needs to be positive or null. When 
no wind is present, the power released by accretion accounts 
for only one third or less of the power available, namely 

Pvis 3 



P acc TZ e ,e V An J \ An J 



(34) 



where the value of lZe,e is set to 3/2 to match the physical 
properties of the external SAD. Thus, energy deposition by 
viscosity (the value of D) is of critical importance only for 
a small radial extenl^] Ari <C n ■ We will come back to this 
crucial issue later. 



2 In standard accretion disc theory, the power lost at the exter- 
nal radius by viscosity is negligible with respect to the accretion 
power Pace because r e 2> n. On the other hand, the "zero torque 
condition" (i.e. T r< f, = 0) assumed at the inner boundary ri for- 
bids any incoming flux of energy. Both assumptions (r e 2> n and 
T r ^ = 0) are violated in the context of fan-shaped winds. 



The luminosity 2P laa of a thin disc is the power that has 
been dissipated into heat by turbulence inside the disc and 
radiated away at its surfaces. In a steady-state framework it 
writes 



2P rs 



dV 



pvv I r— \ + r] m J^ + rj^Jj, 



(35) 



where the integration is made over the volume V of the disc. 
We can rewrite 2P ra( j = f^" Q(r)dr where Q = Q v + Q m , 
namely the viscous term is 



Qv{r) 



GMM a 9 

7 2 W e 



and the magnetic (effective Ohmic heating) is 



, , GMM a 



with T = lZ m E + 



(36) 



(37) 



n ra\_n 



and IZn 



TZeVm is the mag- 



netic Reynolds number (vertically averaged). The factor p.T 
is estimated to be at most of the order of e = h/r in all 
situations. As a consequence, the magnetic dissipation term 
is always negligible with respect to the viscous dissipation. 
Since we are dealing with a thin annulus where most of the 
dissipation is expected to occur at the innermost radius, we 
derive the following estimate for the disc luminosity 



2P rad ~ Q v (n)An 



+ (38) 

z zr° ri 



It is interesting to compare this power to the incoming one 
by viscosity, namely 

2P lad ^ 3 An 1+23 
" 2 n 



p . 

i vzs 



(39) 



In the absence of any ejection (/ = 0), the energy balance 
of the annulus would be simply P acc + P V i S = 2P la d, which 
provides T> — 0. In order for winds to be present (/ 7^ 0), 
the difference factor T> must be positive. 



2.5.4 The power available for winds 

The available power feeding cold fan-shaped winds is 

2Pwznd — Pace + Pvis 2P Taa \ 

GMMg, e fAn , | 3 / p An \ 
2n \ ri ' iz e ,e V n J 

9 An 1 + V ' 



(40) 



2 n Tie, 

where Ohmic dissipation has been neglected. Equating this 



estimate with expression ( 26 1 and using the boundary condi- 
tion 1Z, 



qp 2 f 



e,e — 3/2 provides the following important constraint 
An 



23- 



(41) 



that must be fulfilled for any kind of fan-shaped wind. This 
expression relates the mass loading parameter ft of fan- 
shaped winds to disc physics (q, f) and in particular to the 
difference factor T>. The average magnetic lever arm of fan- 
shaped winds writes J = 1 + 23// while their power is 
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2P„i n d — 7. — (22? — /) 

2n 



(42) 



and is directly controlled by 23. The value of / must be 
obtained by a complete calculation of the two-dimensional 
axisymmetric accretion-ejection problem. In the absence of 
such calculation, we use / as a free parameter, simply ful- 
filling / < 223. It is clear that, unless T> is of order unity or 
larger, fan-shaped winds are just epiphenomena with respect 
to accretion, both in terms of power and mass loss rates. As 
a consequence, this implies that the power feeding such fan- 
shaped winds must be entirely supplied by viscous transport 
from the regions below ri. It is important to remember this 
point. 



3 THE VARIOUS KINDS OF FAN-SHAPED 
WINDS 

In the previous section we investigated some generic prop- 
erties of the disc and its wind, regardless of the configu- 
ration. We just assumed the existence of steady-state cold 
fan-shaped winds launched from a disc annulus of extent 
Ar; ~ h. This annulus is settled between an inner radius 
ri where the poloidal magnetic field is near equipartition 
~ 1) and an outer radius r e , where a standard accre- 
tion disc is established with <C 1 and !Z e ,e = 3/2. In the 
following, we analyze the three possible configurations en- 
visioned in the introduction and discuss also the particular 
case of the X-wind. 



3.1 T-winds 

A Terminal-wind is located at the outer part of a JED, a 
disc which is launching cold jets. By definition, one expects 
p = B^/Bz of order unity and thus TZ m ,i — pr/h ~ e _1 
(see Appendix |Aj. As a consequence, the Reynolds number 
at the JED interface is TZ e ,i = IZmjPm 1 — Pm lr /h, which 
translates into 



23 



3 h 
2 r 



(1 - f)V„ 



(43) 



A turbulent disc with P m ~ 1 would display 23 ~ — 1 and the 
situation considered here would be energetically impossible. 
The reason is quite obvious: the annulus would loose much 
more energy at r e by viscous friction than what it would 
gain at n, even considering the case / = 0. Thus, either the 
inside-out transition from a JED to a SAD is done over a 
larger extension, namely Art ~ n or more (but then the 
assumption of winds taking the shape of a fan would break 
down), or P m is larger than r/h so that 23 > becomes non 
negligible. 

This has however strong implications on the inner JED 
physics. Indeed, a JED requires m™ 9 ~ a m ,i of order unity 
( |Ferreira|[l997| |Zanni et aL||2007[ ), which implies that the 
ratio of the magnetic to the viscous torque becomes 



A, 



(T m e 



(44) 



since C = 3/2+3£ in a JED where M a oc r 5 . Thus, P m ~ r/h 
or larger implies a viscous torque that remains comparable 
to the magnetic torque in the whole JED zone. This is why 



some power can be successfully conveyed outwardly by vis- 
cosity and ultimately feeding the T-winds. Defining ri n as 
the innermost radius of the JED, this power can be esti- 
mated as 



2P W 



^r^\ (21 



2D-f) 



f 



(45) 



where P^cc.jed is the available power for the JED and r» 
corresponds to its outer radius (rj, Fig. 1). If P m ~ e~ is 
such that 23 is of the order of — ^ only, then T-winds would 
be just epiphenomena located outside a more powerful disc 
wind and allowing the transition from an inner JED to an 
outer SAD. This sounds quite reasonable and might deserve 
some attention. But if Vm is larger so that 23 if of order unity 
then the inner JED undergoes a significant loss of power by 
friction at its outer radius and one cannot neglect this loss 
in the global energy budget of JEDs. 

However, regardless of these considerations, Vm ~ £ _1 
or larger would require the following turbulence parameters 
in the disc 



a v ,i ~ e 
a' m ~ 1 

Urn ~ 1 



(46) 



Obtaining such a large viscosity is quite problematic and is 
probably an indication that one or more assumptions are 
actually violated. 



3.2 M-winds below R co 

When the stellar magnetic field is such that R t < n < R co , 
namely a disc truncation below the co-rotation radius, ac- 
cretion onto the star is able to proceed along magnetospheric 
curtains (accretion columns or funnels). This is possible be- 
cause the stellar magnetic field brakes down the disc mate- 
rial below r^, in a way that is pretty much the same as in 
a JED. While the stellar magnetosphere is strong enough 
to truncate the disc right below n with fit of order unity, 
it efficiently brakes down the disc material giving birth to 
a sonic accretion flow at ri so that an accretion funnel flow 
can take place ( Romanova et al.|2002 2009; Bcssolaz et al 



2008 Zanni & Ferreira 20091. We identify therefore n as 



the base of the funnel flow (termed rbf in Bessolaz et al. 
2008 1 . This situation is therefore very similar to the T-wind 



case above, where the two torques (viscous and magnetic) 
are braking down the disc plasma. 

If cold winds are to be launched between n and r e , then 
the poloidal magnetic field must be bent enough so that 
7Z m ,i ~ r/h, leading to a Reynolds number 1Z e ,i ~ P^r/h 
as in the case of T-winds. The conditions for these M-winds 
are therefore exactly the same as for T-winds, namely P m > 
r/h. Since m s ,i ~ 1 here also, the turbulence parameters 



must also verify ( 46 1 



3.3 M-winds beyond R co 

When the stellar magnetic field is such that R co < R t < 
ri, then no accretion columns can be formed. Indeed, the 
magnetosphere below n is actually accelerating azimuthaly 
the disc material, which is pushed outwardly. In other words, 
the magnetic stresses due to the closed magnetospheric field 
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lines deposit stellar angular momentum into the disc and 
accretion cannot proceed anymore below n. A steady-state 
configuration can be maintained only under two additional 
conditions: 

(i) Mass conservation requires the complete deviation of 
disc material into the M- winds, namely / = 1. 

(ii) Viscous stresses must be able to transport the stellar 
angular momentum radially beyond r\ . This translates into a 
viscous torque that must be positive between rt and r e . If it 
were not for the torque due to the M-wind, accretion would 
stop. Assuming that the disc recovers a standard structure 
at r e demands a transition zone where the negative torque 
due to the M-wind is larger than the positive torque due to 
the viscosity. 

Since there is no mass accreting below rj, the power 
released by accretion is negative 



p - 



GMM ae (Ar. 



GMM ae 

2n 



(47) 



(48) 



2r; \ r 
and one gets simply 

2P W ^^^(2Z>-1) 

or, equivalently, winds with an average magnetic lever arm 
J — 1 + T> with T> > 1/2. The power feeding the jets is here 
provided by the viscosity but its source is the rotational 
energy of the star. 

The difference factor can also be written 

23 = |i2M-l (49) 

which is always valid, including in the case / = 1 considered 
here. The presence of mass Ei stuck at Ti (with a vanishing 
radial velocity) is required for the viscous transfer of energy 
from the rotating magnetosphere below Ti to the annulus of 
extent Ar^. M- winds of this kind could therefore only exist 
if 

" (50) 



> 



2 Ei 



Clearly, Ei cannot be much larger than E e otherwise the 
thin disc (near Keplerian) approximation would break down 
as this would generate a significant radial pressure gradient. 
It is therefore safe to argue that the viscosity coefficient a v ,% 
at Ti must be at least comparable to the one at r e . 

Mass flux conservation requires 2p + u+ Ari ~ E e u e in 
a steady state framework, where all accreted mass is being 
deflected into the winds. Given our geometry with An ~ 
h, this implies u e ~ uf ~ C s in order to start a wind. 
In other words, the disc plasma is entirely deflected in the 
vertical direction into a wind only if the incoming velocity 
is already close to the sonic speed. If this were not the case, 
mass would simply accumulate at the magnetopause Ti and 
the steady-state approximation would break down. Because 
u e — mD^Cs = 3a ViE eC s /2, the condition u e ~ C s cannot 
be fulfilled unless a„, e ~ e~ . 

In order to halt accretion at ri (m 3] i ~ 0), the viscous 
torque must balance the magnetic torque due to the M- 
winds, namely be such that ~ — 2qifii/a V: ie < 0. Using 
Eq. ( 15 1 in the transition zone with f — 1 gives 



An 



Tie 

7^ 



- 2 



-3 



3a v 



Ari \ 2m. 



- 2 



(51) 



The fact that the Reynolds number decreases drastically 
when approaching r» is consistent with a viscous torque be- 
coming positive. The sonic Mach number consistent with the 
disc angular momentum conservation writes 

(52) 



22+| 



h 



in a region where Ar; ~ h and fii ~ 1. If a Vl i > qi then 
m s ,i ~ a v ,iE > a„, e £ ~ 1- This is inconsistent with an ar- 
rested accretion. The only alternative is therefore to have 
a v ,i <C qt which provides m 8 ,» ~ a^ i q~ 1 e that can be in- 
deed vanishingly small. However, this requires qi 3> e _1 , 
namely magnetic field lines dramatically twisted. Note that 
it is because the viscous torque is huge at r e that the toroidal 
magnetic field must be so large. 

The turbulence parameters must verify in this case 



(53) 



which is again troublesome, regardless of the value of Vn 



3.4 M-winds from R co '- X-winds 

The limiting case when R t < r, ~ R co is interesting be- 
cause it may possess the good sides of the above two M-wind 
types without their drawbacks. Such a situation is commonly 
known in the literature as X-winds. First, some rotational 
energy might flow from the star to the disc beyond rt thanks 
to the viscosity (see discussion p786, left column, last para- 



graph in Shu et al. 1994a I. Second, the ejected fraction / 
could remain smaller than unity because the disc below n 
would have the ability to form accretion columns. Let us 
then investigate this situation. 

The expressions given in Sect. ( |2.5.4[ ) are valid in the 
case of X-winds as well and can be cast into the following 
form: 



J — 1 _ Ti V 

~ q Ar~ ~ ]Pf 



(54) 



with the constraint / < 2T>. So far, / (or T>) is unspecified 
and free as long as a full MHD calculation including the disc 
remains undone. However, actual explicit calculations of X- 



winds do exist and provide useful informations ( |Shu et al. 
1994b||Najita fc Shu|1994||Shu et al.|1995||Shang et al.|1998| 



2002 



Cai et al.||2008 1. These computed 2D flow structures 



are fan-shaped cold winds expanding from a point source 
with an initially super slow-magnetosonic speed, matched 
to an asymptotic jet cylindrical solution. 

All published solutions display j3 of order unity (be- 
tween 1 and 3) and J between 2 and 7 (see e.g. Table 4 in 
Najita fc Sriu]|1994| and Table 3 in |Cai et al.1|2008| |. As the 
underlying disc is a mere boundary condition, some assump- 
tion has been made in order to derive the value of the ejected 
mass fraction /. Indeed, it has been assumed that X-winds 
transport away the exact amount of disc angular momentum 
accreted onto the star, namely M a jf2jrf = 2M w JQjr 2 . This 
translates into fj = 1 - / (Eq. (4.7a) in |Shu et al~1|l994a 
gives a more complete relation, but published models satisfy 
this last scaling) or O = 1 — 2/ with / ranging from 0.1 to 
0.4 (Table 5 in|Najita fc Shu||1994[). Thus, X-winds as they 
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appear in the literature require 

An 



V = 1 - 2/ ~ 1 and q. 



(55) 



Note that these values are consistent with those found in 



the aforementioned papers: see for instance Eq.(4.19) in Shu 
et al. ( 1994a[ ) providing the power carried by the wind and 
the value assumed for q in |Najita fc Shu ( 1994 1 (see also 
discussion §2.6 in |Shu et al.||1994b |. 
Now, Eq.( 15 1 writes for any / < 1 



= 



v 



An V 1 - / 



-/ 



An 



(56) 



which shows that the viscous torque is indeed positive as in 
M- winds beyond R co . From the definition of the difference 
factor T> one gets using 7£ e , e = 3/2 



1Z e 



3 1-/ _ m s ,i 
2 1 + V ~ a v ,i£ 



(57) 



showing that the sonic Mach number at r» must be m 8 ,i = 
2gi/Xi + (,iOt v ,iE = 3a v>i e/4: for T> = 1 — 2/. With the con- 
straints that both [ii and T> being of order unity and q ~ e, 
one gets 



(58) 



where the last condition has been derived using lZ m ~ 
e -1 (cold ejection). These transport coefficients are those 
required within the disc so that X-winds can indeed be 
launched with T> — 1 — 2/ ~ 1. The almost perfect match- 
ing between the two modes of angular momentum transport 
(turbulent /radial, wind/vertical) gives rise to a slight unbal- 
ance, allowing henceforth accretion. 



4 DISCUSSION 

We have shown that, under reasonable circumstances, all 
fan-shaped winds (T-winds and all kinds of M-winds, in- 
cluding X-winds) put severe constraints on the underlying 
MHD turbulence. It appears that our conventional view of 
turbulence can be hardly reconciled with any kind of fan- 
shaped winds from thin (Keplerian) discs. Before discussing 
further the consequences of our findings, let us first address 
some caveats of our analysis. 



4.1 Caveats of our analysis 

4-1.1 MHD turbulence and anomalous transport 

In order to design steady-state models of accretion discs, one 
needs to take into account the turbulent transport of both 
angular momentum and magnetic field. It has been assumed 
that turbulence leads to a local transport that has been in- 
cluded in our equations as an anomalous viscous torque and 
an anomalous resistivity term. Self-sustained turbulence in 
accretion discs is thought to be triggered by the magneto- 
rotational instability (Balbus & Hawley 19911. Thorough 



instability does indeed provide an outwardly directed angu- 
lar momentum transport that, indeed, can be crudely de- 
scribed by a local turbulent viscosity in a vertically aver- 
aged disc v v (see Balbus 2003 ) . Many important results have 
been obtained in the context of shearing boxes (e.g. Pessah 
et al.|2007 Lesur fc Longaretti|2007 and references therein). 
One of these results concerns the existence of an anomalous 



magnetic field transport in MRI-driven turbulence ( Guan & 



Gammie||2009l |Lesur fc Longarettil[2009l |Fromang fc Stone 



2009 ) . While previous works were only focused on analyzing 



numerical and analytical analysis tend to show that this 



the angular momentum transport, these studies investigated 
the magnetic field dissipation. They concluded that it could 
also be roughly described by some local Ohmic resistivity 
term with an effective magnetic Prandtl number V m of or- 
der unity (Lesur & Longarctti 20091. 

Two words of caution though. First, these studies were 
done only with very weak fields, namely [t, <C 1 and it 
remains therefore to prove that the behavior and scalings 
found would remain valid for a magnetic field near equipar- 
tition. Second, it is not clear how the boundary conditions 
in shearing boxes affect the development of the turbulence. 
Fully resolved and converged 3D simulations are therefore 
needed to verify if some non local behavior is present or 
not (see e.g. |Beckwith et al.|20 09, 2011}). Anyway, given our 
present knowledge, we believe that using the functional form 
of the viscous stress tensor for the turbulent transport with 
a constant V m throughout the disc remains a reasonable, 
albeit limited, approach. 



4-1-2 Cold wind approximation 

In our analysis, we neglected the enthalpy in the disc en- 
ergy budget and assumed field lines inclined by more than 
30° with respect to the vertical. Both assumptions are not 
crucial but allow a simpler treatment of the equations. The 
amount of enthalpy carried in by the flow (both in the jet 
and at rj) remains negligible as long as the disc is geometri- 
cally thin. The thin disc approximation goes together with 
the near-Keplerian assumption and provides indeed a lim- 
itation to our analysis. This work must then be seen as a 
first attempt to understand whether fan-shaped winds can 
be steadily launched from near Keplerian discs. 

In fact, including thermal effects in the launching pro- 
cess should not change our findings. Indeed, what is critical 
here is the strength of the poloidal magnetic field and its 
geometry. The magnetic field strength is assumed to decay 
outwardly so that /x goes from about unity at n to a tiny 
value at r e . Correspondingly with this decay, the poloidal 
field experiences a quasi-spherical dilution so that ~ B z 
between n and r e is a natural outcome due to the geometry. 
Maintaining this geometry despite the natural tendency of 
the field to diffuse outwardly is an open issue that we do 
not address here. We assume a fan-like geometry and look 
at the dynamical consequences on the disc. 

One could argue however that, contrary to the T-wind 
case, M-winds are located just outside a closed magneto- 
sphere. This means that below r% the magnetic field lines are 
actually directed towards the star. Thus, must undergo 
a steep transition from ~ / Hoh ~ B z /fi h > be- 
tween n and r e and a negative value below n. The poloidal 
field lines must then be almost vertical somewhere around 
n, defining thereby a zone from where no cold ejection could 
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take place (a "dead zone", Shu et aT]|1994a |. However, this 
will affect only few field lines around the radius where 
vanishes and r, must then be seen here as the radius where 
the fan-like geometry has already provided the correct bend- 
ing. 

4.2 Are stationary fan-shaped winds possible? 

4-2.1 T-winds and M-winds below co-rotation 

Both situations are very similar and energetically forbidden 
for Vm ~ 1. On the other hand, even the case V m ~ e _1 is 
problematic as it would require a huge viscosity coefficient, 
ct v ,i ~ e . This is clearly against our current knowledge 
of MHD turbulence, where a v < 1. Nevertheless, it may 
be of some interest to make numerical experiments, with 
a prescriptions and these parameters, to see whether fan- 
shaped winds would be launched. 

Our guess is however that our assumption of the very 
existence of fan-shaped winds is violated. Precisely, that no 
JED-SAD or Magnetosphere-SAD transition giving rise to 
winds can be done in a radial extent Art ~ h. But since 
transitions like those are anyway necessary, we conjecture 
that they must be done on a larger scale Ar^ ~ r, for a 
steady-state to settle. 



4-2.2 M-winds beyond co-rotation 

This configuration is even more problematic. It not only 
requires a v ~ e _1 as the other two, but a steady ejection of 
all the incoming mass requires a huge magnetic shear q 3> 
e - . It is doubtful that such a large winding of the magnetic 
field lines could be steadily maintained. Reconnection will 
probably take place so as to (violently) relax the magnetic 
field configuration to a lower energy state. Alternatively, for 
instance if a v < 1, then the disc material accreting below r e 
would be unable to be vertically deflected with / = 1 and 
mass would start to accumulate. 

Our guess is that steady state M-winds will never be 
established beyond R co - 



where i/ n i is the neutral-ion collision frequency and p ~ 
Pn + pi = p«(X + X) with p n (pi) the neutral (ion) density 



and X the ionization fraction ( Ferreira fc Pelletier||1993a l 
Normalizing this diffusivity by Qnh gives another alpha 
coefficient, measuring the level of the ambipolar diffusion 



CtAD 



(1 + A) 



(60) 



with T ni — v~l ■ It is clear that a fluid approach requires 
^■KTni <C 1 within the disc and it is indeed the assump- 
tion made in X-wind theory. However, this implies ctAD ~ 
pflxTni <C /i < 1, much smaller than the value a' m required 
to maintain q ~ e. We thus believe that, as long as some 
MHD turbulence is triggered and sustained within the inner 
disc, it will be the dominant diffusion term (note also that in 
the ambipolar diffusion regime, no MRI would be permitted, 
Bai fc Stone|2011| >. 

A second difference is in the use of the disc angular mo- 
mentum. While we explicitly used it to compute the sonic 



Shu et al. (1994a I, namely m 



Mach number m s (Sect 2.31, we found only an order of esti- 

Qi<T n i <S 1. How- 



mate in 

ever, the actual value of m s remains unspecified. Within our 
framework and with 23 of order unity, we obtain m s ,; ~ e 2 , 
which is indeed accretion done at a very slow pace. This is 
understandable as the angular momentum removal by the 
wind has to overcome the viscous deposit. But this is how- 
ever problematic as the formation of accretion columns does 
require an almost sonic accretion at the base of the funnel 



flow (Bessolaz et al. 20081. How will the disc material be 



loaded onto the magnetospheric field lines below r; remains 
an unsolved issue within the X-wind scenario. 

Another drawback concerns the value of the turbulence 
parameters that are required for powering X- winds with 23 ~ 
1 in a steady state. The low accretion speed combined with 



the small magnetic shear require (see Eq 58 1 



(61) 



namely a huge anisotropy in magnetic field dissipation, re- 
gardless of the value of Vm- This is inconsistent with our 
current knowledge of turbulence in shearing boxes. 



4.2.3 X-winds 

From the above discussion it appears that the only remain- 
ing configuration that could lead to stationary fan-shaped 
winds is when is near the co-rotation radius, namely X- 
winds. As already stated, we recovered with our formalism 
the main results of X-winds concerning the strength of the 
magnetic field (/i; ~ 1, §2.8 in |Shu et~a l. 1994a), the value 
of the magnetic twisting q ~ e, the link between the MHD 
invariants (/?, J) and the global energy budget (23, /). There 
are however two main differences. 



First, Shu et al. (1994a) do not use the usual single 



fluid MHD induction equation for the magnetic field. In- 
stead, ambipolar diffusion is assumed to be the main source 
of magnetic diffusion, with a possible comparable role of the 
Ohmic resistivity (see their §3.1). However, the ambipolar 
diffusion term can be translated in order of magnitude into 
a diffusivity 



VAD 



V li 1 



(59) 



4-2.4 Designing different X-wind solutions 

The double goal sought by the X-wind model, as theorized 
in Shu et al. ( 1994a I, is to explain YSO jets while simultane- 



ously braking down the protostar. More precisely, X-winds 
should be the dominant jet (in terms of appearance, mass 
flux and power) while carrying away enough angular mo- 
mentum so that the star is not being spun up by the accret- 
ing material (zero net torque condition). These requirements 
lead to the choices of 23 = 1 — 2/ ~ 1, q ~ e, J of order of 
a few and p of order unity. However, these values raise the 
issues discussed above. 

The difficulties faced by these M-winds can be under- 
stood in a broader framework. As expressed by Eq. ( 33 1 23 



measures the net energy deposition by viscosity in the zone 
of extent At\ ~ h. This deposit will eventually power the 



jets (Eq 42 1 and that is why M-winds of any kind require 23 



of order unity to be of any significance. The problem is the 
tiny width of the launching zone, which imposes that 23 ~ 1 
requires a Reynolds number of order unity as well. Other- 
wise there is much more power lost at r e than the one gained 
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at ri and T> ~ — 1, Eq.(57l. This translates unavoidably into 
m s ,i ~ Qu.iE which is of order unity only for a v ^ ~ e -1 . 

Let us design new M-wind solutions by demanding 23 ~ 
1 and m Sji ~ 1 only. 23 of order unity requires fj ~ — and 
the viscous torque is accelerating the disc material. Then, to 
maintain accretion at a sonic pace, the magnetic torque must 
be correspondingly huge with q± ~ -^ r . As a consequence, 

— - it one 



541 imposes (3 



the interplay with the jet (Eq I 
desires winds with asymptotic speeds comparable to those 
observed. In practice, this means that new super-fast MHD 
wind solutions must be looked for. 

Nevertheless, even if such wind solutions could be 
found, they will demand the following a coefficients in the 
disc 



a 



V.l 



(62) 



with both V m ~ e~ and \ m ~ e~ . This is again trou- 
blesome, not only because of the strange anisotropy that 
would be required but because of the overwhelming mag- 
netic compression on the disc implied by q ~ e _1 with a 
near equipartition field: no vertical balance could be found 
(see Appendix [B|) . 

Finding M-wind solutions that would be consistent with 
our conventional knowledge on MHD turbulence, namely 
Vm ~ Xm ~ 1 is impossible. Indeed, 23 > (even small) re- 
quires TZ e ,i of order unity while the geometry itself requires 
TLm,j ~ e~ x . Thu s, the minimum requirement is Vm ~ e" 1 
( Shu et al.||2007[ ), namely a m ,i ~ 1 and a Vl i ~ e . Now, if 
one req uires that q is at most of order unity with 23 ~ 1, then 
Eq.(54l requires f/3 2 ~ e. Thus, new dynamical flow calcula- 
tions should be looked for since, for J ~ 1+23// comparable 
to values used in current X-wind models, the mass loading 
parameter should be e 1//2 times smaller. How this impacts 
on the wind dynamics needs to be investigated as j3 <C 1 
describes a priori a matter dominated flow (Rbp ~ 1) and it 
is not clear that any super- Alfvenic wind could be obtained. 



4.3 Comparison with numerical simulations 

4.3.1 T -winds 

There is no proper ab initio MHD simulation of turbulent ac- 
cretion discs offering a radial transition from an inner JED to 
an outer SAD. Simulations with well defined JEDs are those 



done using a prescriptions (Casse & Keppens 2002 Zanni 
|et al.|2007|[Tzeferacos et al.|2009||Murphy et al.|201o| . How- 
ever, in most simulations \x is a constant of the radius and 
the outer radius of the JED is continuously increasing in 
time. The only simulation with a steady-state super-fast jet 
launched from a zone of constant extent has been performed 



by Murphy et al. (20101, with an initial magnetization de- 
creasing with the radius. Clearly, there is no fan-shaped wind 
making the transition from the inner zone with super-fast 
jets to the outer zone embedded in a current-free magneto- 
spheric corona (see e.g. Figs 1 and 2 in Murphy c t al.|2010 |. 
It might be because V m is of order unity in these simula- 
tions: in that case, no sharp transition (i.e. Ar^ ~ h) with 
a spherical dilution can be energetically feasible. Instead, a 
zone of sub-fast but still super- Alfvenic jets is surrounding 



the inner zone from n ~ 5 to r e — 13, namely with an 
extent Ar^ ~ r». Alternatively, it may be just an effect of 
the initial magnetic field distribution, since the simulation 
lasted for 953 inner Keplerian orbits only, which is roughly 
the accretion time scale from r = 13. Performing longer sim- 
ulations is a task that should deserve some attention. 



4.3.2 M-winds 

MHD simulations of a star-disc interaction with an axisym- 
metric dipole field have been done by several authors, mostly 
using a prescriptions for the viscosity and magnetic diffu- 



sivities (see e.g. Romanova et al.|2002|[2009 Bessolaz et al 



2008 Zanni & Ferreira 2009 and references therein). How- 



ever, all of them were done with Xm = 1 and with both 
a v and a m smaller than unity. As a consequence, they are 
all outside the parameter range required to produce steady 
state M-winds as described in this pape r, even if some of 
them were performed with V m ~ £ _1 (Zanni & Ferreira 



2009 Romanova et al. 


2009). 


Zanni & Ferreira 


(2009) 



( 2009 1 showed that decreasing the 



value of a m = a m gives rise to an opening of the stellar mag- 
netic field and the loss of any causal link between the star 
and the disc. This is because diffusion is not efficient enough 
and the star-disc differential rotation leads to a huge gener- 
ation of toroidal field. The response of the magnetosphere to 
this increasing shear is to inflate and reconne ct (|Romanova| 
et al.||1998| and discussion in |Matt fc P udritz 200_5j), break- 
ing thereby the connection with the star. Depending on the 
parameters, some unsteady ejection events (plasmoids) can 
be both thermally and magnetically launched after a recon- 
nection. But these plasmoids carry no self-confining electric 
current and have nothing in common with steady state M- 
winds. The existence of steady "conical winds" was how- 



ever reported in the simulations of Romanova et al. ( 2009 \ 
when the disc is truncated below the co-rotation radius. Such 
winds seem to require V m ^ 1 to exist but the mass loss / 
decreases when Vm increases. Although the origin of these 
winds has not been fully addressed, they clearly do not re- 
semble fan-shaped winds. Instead, they look better like usual 
disc winds that would be launched only from a small radial 
extent (the region where fj, < 1) and confined by the pres- 
ence of an outer open magnetic flux. When the disc is trun- 
cated beyond the co-rotation radius, the system enters the 
so-called "propeller regime" (Romanova et al. 2003 2005 
Ustyugova et al.||2006 and references therein), which is un- 
steady. 

Ab initio simulations that compute also the disc turbu- 



lence have been performed recently, both in 2D ( Romanova 
et al.||2011[ ) and 3D ( |Romanova et al.||2012[ >. These simula- 



tions are highly valuable as they are not biased by any a 
prescription. A first important result is that the disc trun- 
cation radius is indeed located at the position predicted us- 
ing analytical arguments, namely where /t ~ m s ~ 1 (see 
Bessolaz et al.|2 008). A second result is that, in the 2D sim- 
ulations where the effective magnetic resistivity seems to be 
very small (V m S> 1), still no M-winds were obtained. This 
holds true in 3D ( Romanova et al.|2012 \ . This deserves fur- 
ther investigation, but the reason might be that a turbulent 
disc structure is unable to provide the required level of tur- 



bulent transport as measured by the conditions ( 62 1 . 
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5 CONCLUSIONS 

We addressed in this paper the existence of cold, steady- 
state, fan-shaped winds arising from near Keplerian discs. 
Such winds are assumed to be launched from a disc annulus 
of radial extent Ar of the order of the disc thickness 
h. They may play the role of a buffer between different 
magnetic regions in protostellar accretion discs. There are 
two main configurations (see Fig. 1). The first one, referred 
to as " Terminal- winds" (T- winds), acts as a transition 
layer between an inner Jet Emitting Disc and an outer 
standard accretion disc. The second kind of fan-shaped 
winds are the "Magnetospheric- winds" (M- winds), linking 
an inner closed stellar magnetosphere to the outer standard 
accretion disc. Although described by the same set of 
MHD equations as usual disc winds, fan-shaped winds have 
interesting properties (such as the quasi-spherical dilution 
of the magnetic flux) that make them worth to be analyzed. 



The constraints imposed by all published super-Alfvenic 
X-wind flow models on the underlying disc structure lead to 
quite an extreme transport coefficient set-up. For instance, 
the degree of anisotropy of the magnetic field dissipation 
must scale as v m /v' m ~ (h/r) 3 . 

None of the various transport coefficients set-ups 
obtained for T and M-unnds are in agreement with our 
current knowledge of MHD turbulence in near Keplerian 
discs. One key assumption made in this work is the small 
radial thickness Ar ~ h of the wind launching zone. But 
if Ar ~ r, the geometry would not be consistent anymore 
with a spherical dilution of the magnetic flux, which is the 
cornerstone of fan-shaped winds. Thus, unless MHD tur- 
bulence at a JED-SAD and Magnetosphere-SAD interface 
has quite peculiar properties, we believe that steady-state 
T-winds and M-winds of any kind (including X-winds) will 
not be realized in Nature. 



Taking into account the physical conditions imposed 
by both inner and outer boundaries and considering the 
whole set of resistive, viscous MHD equations, we were 
able to describe for the first time the mass, energy and 
angular momentum balance of the near Keplerian disc. To 
drive winds with a significant power, the annulus must be 
fed at its inner edge with more energy than what is lost 
from winds, radiation and viscosity at its outer edge. We 
found that the only possible way to power such winds is to 
allow for an efficient viscous energy flux at the inner edge. 
Such influx of energy by viscosity goes along with an influx 
of angular momentum. Thus, for accretion to take place 
despite this angular momentum excess, the torque due 
to the fan-shaped winds must be dominant wrt the local 
viscous stresses. Now, both the magnetic geometry and 
the balance of these angular momentum/energy fluxes put 
stringent constraints on the underlying MHD turbulence in 
the disc. 

These conditions were analyzed using local a prescrip- 
tions for the transport coefficients: viscosity i/ v , magnetic 
diffusivities u m and v' m . We found the following results: 

(f) For the case of T-winds and M-winds whose inner 
radius r< is located below the co-rotation radius R co , we 
found that fan-shaped winds can exist only if the effective 
magnetic Prandtl number is larger than n/hi (hi being the 
disc thickness at r») and a viscosity v Vi i larger than flK,irihi 
(a v ,i > n/hi). 

(2) In the case of M-winds with R co > n, we found a 
similar viscosity amplitude and a magnetic field configura- 
tion likely to be non-steady, as it exhibits a huge magnetic 
twisting prone to reconnection. 

(3) We also considered in our study the particular 
case of M-winds whose internal radius coincides with R co , 



namely X-winds as first proposed by Shu et al. (1994a I. 
We argue that, contrary to the assumption made in the 
earlier X-wind papers, ambipolar diffusion is probably not 
the dominant source of magnetic field diffusion at the star- 
disc magnetopause. Instead, since the large scale vertical 
magnetic field that is threading the disc must be close 
to equipartition, we expect strong field instabilities to be 
triggered there and to saturate into a self-sustained MHD 
turbulence (described here as local transport coefficients). 



While stationary T-winds and M-winds appear to be 
ruled out, would such type of outflows still be possible if 
the steady-state assumption is released ? Let us for instance 
consider that a thermal-viscous instability is triggered in 
the disc (see Lasota|2001 for a thorough review on the Disc 
Instability Model). Its outcome is a thermal front, of width 
Ar ~ h, which propagates at almost sonic speeds, namely 
Ufront ~ ctvCs, where a v could be as high as 0.1 (Meyer 



1984 Menou et al.|1999 l. Such a large accretion speed could 
advect large scale magnetic fields in the disc and give rise 
to a magnetic configuration as envisioned for T-winds. One 
could then wonder if fan-shaped winds just require specific 
conditions inside the disc, such as those necessary to trigger 
this instability. But under these circumstances, not only 
stationarity but also several other of our assumptions break 
down. Indeed, steep radial gradients of disc column density 
and temperature, both associated with this thermal front, 
can potentially change the rotation and specific angular 
momentum profiles so as to allow for a local Rayleigh 



instability (Lin et al. 19851. While our present analysis of 



fan-shaped winds took advantage of the existence of MHD 
invariants (and the strong links they provide between the 
winds and the disc), none of these would exist anymore in 
an outbursting disc prone to a thermal-viscous instability. 
Such links would disappear as the time scales involved 
(front propagation, wind readjustment) become too close. 
Note also that a magnetic configuration suited for M-winds 
could in principle be obtained in a time-dependent way. 
Indeed, it could be achieved if a large scale magnetic field 
advected towards the central object is suddenly halted by 
the magnetic pressure of the closed stellar magnetosphere 
and bounces back. These unsteady situations are however 
beyond the scope of our analysis, as full 2D time-dependent 
simulations would be needed to investigate such dynamical 
situations. 



Actually, unsteady ejecta are indeed observed in many 
MHD simulations of star-disc interaction. Both in those 
computing the disc (Romano va et al.||2009 20111 or where 
it is only a platforrnfFendt fc Elstner||2000| |Fendt||2009 1. 
As argued in Ferreira et al. ( 2006[ ), such ejecta are actually 
welcome. Indeed, they could produce time-dependent events 
that would be confined and channeled by an outer steady 
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disc wind (driven from a JED), possibly explaining both jet 
variability and the presence of knots. Also, these ejecta could 
transport a sizable fraction of the stellar angular momentum 
and provide thereby a natural way to brake down the cen- 
tral protostar (see e.g. Ferreira et al.|2000 Zanni & Ferreira 
|2009[ ). According to this picture, there will be no unique 
MHD model capable of explaining all observational features 
of YSO jets. These would be instead a multi-component flow, 
made of a central core filled in by the stellar wind, an inter- 
face layer of unsteady plasmoids and an outer disc wind as 
the ultimate confining agent. 
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APPENDIX A: THE POLOIDAL MAGNETIC 
FIELD AT THE DISC SURFACE 

An axisymmetric magnetic bipolar topology is described by 
a poloidal field B p = ^Va x e^,, where a(r, z) is an even 
function of z, and an odd toroidal field B^. A magnetic sur- 
face anchored at a radius r is a surface of constant poloidal 
magnetic flux that can be labelled by a(r, z) — a (r , 0). The 
bending of the field lines at the disc surface (z ~ h) can be 
evaluated as 



h 



dz 2 



( da \ ' 



fda \ ' 



(Al) 



where I is the vertical scale of variation of the magnetic flux. 
The radial distribution of magnetic flux through the disc is 
an unknown function. However, by definition, fan-shaped 
winds assume a steep decrease of B z with the radius, on 
a scale Ar ~ h. This is in strong contrast with extended 



disc winds as 



Ar ~ h <C r. As a consequence 



da /dr 

B^/B z ~ hAr/l 2 ~ h 2 /l 2 . Extended disc wind configura- 
tion would lead to a vertical sc ale variation of the magneti c 
flux such that B+/B z ~ hr/l 2 (|Ferreira fe Pelletier|l993a|. 



Cold fan-shaped winds are launched by the very same 
magneto-centrifugal mechanism than large scale disc winds. 
They must thus fulfill the |Blandford fe Payne| ( [19821 ) cri- 
terion. This geometrical criterion translates into B^ / B z of 
the order unity, which then implies I ~ h. Contrary to usual 
disc winds, B z varies on the disc scale height. However, this 
is not crucial for our estimates since — B z (z — h) is 
smaller than but of the order of B z (z = 0). 

The cold, fan-like geometry of the magnetic field re- 
quires a toroidal electric current density at n 



J<t> = 



dB z 
dr 



B z 



(A2) 



j_ r&B^ 

p a \dz dr ) ~ fj, a h 

where ci > is of order unity. According to Ohm's law pL 
the balance between field advection and diffusion produces 



Jo 



ViUy B z 



71, 



B z 



(A3) 



where lZ m is the effective magnetic Reynolds number. Thus, 
cold fan-shaped winds are possible only if lZ m — ri/h. This 



requires a quite significant accretion speed or, equivalently, 
a sonic Mach number m s = a m 1Z m h/r ~ a m at n. 



APPENDIX B: A NEAR KEPLERIAN 
QUASI-STATIC DISC STRUCTURE 

The steady-state disc vertical equilibrium corresponds to 
a balance between three forces: gravity, magnetic pinch- 
ing force and thermal pressure gradient. Integrating the 
z-component of the momentum conservation equation |2| 
leads to 



Po 



B+ 2 + B 
2p a 

h 



+ 2 



C2 



+ 



pzQ? K 



2/x 
2 dp 



' dz + 



(Bl) 



where C2 > is of order unity. Given the magnetic geometry, 
the magnetic Lorentz force is always negative, so that the 
only force allowing a vertical acceleration within the disc 
is the plasma pressure gradient. This is possible only if the 
magnetic pinching and compression are not too high, namely 
if 



Po 



Po 



1 > 



+ 



<p_ 

B z 



(p 2 + q 2 ) (B2) 



where p = B^ /B z ^ 1 and q = —B^/B z . This translates 
into an upper limit for the disc magnetization, namely p < 2 
in order to satisfy the vertical equilibrium of the disc. The 
bending p of the magnetic field is directly related to the 
magnetic Reynolds number through the resistive induction 
equation p ~ lZ m £ (see Appendix [A]) . A lower bound can 
then be obtained for the magnetic Reynolds number, namely 



Tim > e 



Note that a steady-state situation is achieved 



only if P <g P , ie. p <p„. This simple argument shows 
why a M-wind in the propeller regime cannot be described 
within a steady-state framework. 

In the zone of extent Art, both the bending (B^ /B z ) 
and the shear (B^/B z ) of the magnetic field lines increase 
with distance. But this can be easily compensated with the 
steep decrease in B z (hence p) so that the above condition 
can remain fulfilled. Note that, because of the vertical mag- 
netic force, the real disc scale height is actually smaller than 
the hydrostatic scale height h, as defined in Eq.Q. However, 
this extra compression is at most of the same order than the 
gravitational one so that using h and C s ~ Qxh is good 
enough for our purpose ( Ferreira & Pelletier 1995 Shu et al 
2008b. 



The steady-state angular velocity fi of the plasma is 
given in the disc by the radial component of Eq. Q. Ne- 
glecting the advection term it writes on the disc midplanc 



in 



i + 



dP/dr 
p £l 2 K r 



J<j,B z 
p o 0, 2 K ' 



Cip- 



(B3) 



where C3 is at most of order unity (both signs are allowed). 
The strong radial decrease of the vertical field gives a mag- 
netic support against gravity. However, for the conditions 
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envisioned here, namely fi at most of order unity, a Keple- 
rian rotation law remains a good approximation in a thin 
disc. 



APPENDIX C: THE TOROIDAL MAGNETIC 
FIELD AT THE DISC SURFACE 

The disc induction equation (|4|) for the toro idal component 
can be written as ( Ferreira fe Pelletier|1995 Ferreira|1997 1 



V' m Jr(z) = V' m jr\,-n + ' 



B<hU z 



(CI) 



in the resistive MHD regime, where the last term (advec- 
tion) can be neglected as the vertical velocity is very small 
compared to the toroidal velocity in a thin disc. At the 
disc surface the product r)' m J r must vanish for two rea- 
sons. First, turbulence must decrease so that the flow be- 
comes frozen in the magnetic field. Second, J r itself must 
also decrease to allow for a change of sign of the magnetic 
(wind) torque F^ = J z B r — J r B z . Indeed, from negative at 
the disc midplane, the torque must become positive at the 
disc surface. This is the cornerstone of magnetic launching 
from thin discs JFerreira & Pelletier||1993bl 119951). In the 



absence of dynamo, the only possibility is to balance the 
term rj' m J r at z = with the counter electromotive force 
(cemf) due to the disc differential rotation rB p • Vf2 inte- 
grated over the disc thickness (see for instance Fig. 2 in 
Ferreira fe Pelletier|1993b I . Given the fan-shaped geometry, 
the main contribution to this cemf is the radial one, which 
writes r J B r ^dz ~ — J B r dz ~ —c^licBgh, where 
C4 > is of order unity for a cold wind. Thus, a vanishing 
,l r at the disc surface is realized if, at the disc midplane, it 
amounts 



Jr,o — C4 



Hoh v' m 



C4 B z 



Defining precisely the magnetic shear parameter q as 
B z 1 dB^ 



Jr,o 



I 7 



fi dz 



B+ 

>p_ 

jx h 



(C2) 



(C3) 



we obtain an estimate of the toroidal field at the disc surface 
1 - ^ (04) 

This is an important constraint, known in the field of wind 
launching from thin accretion discs and worth to comment. 
The toroidal field at the disc surface is critical for at least two 
reasons. First, it is responsible for the torque due to the wind 
and determines therefore the disc accretion rate. Second, 
it conveys the power (as the MHD Poynting flux leaving 
the disc) that drives the wind. Its magnitude is therefore 
a key ingredient in magnetic jet acceleration models. The 
above result, due to the induction equation Q, tells us that 
the twisting of the magnetic field is stronger as the MHD 
turbulent diffusivity gets smaller. However, no steady-state 
model can be achieved with g > 1. As a' m decreases, the 
magnetic compression of the disc becomes overwhelmingly 
large (Appendix [B| and the whole structure gets unsteady 
( Zanni et al.|2007[ ). 
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